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Abstract: Apoptotic cell death is coordinated through two distinct (type 1 and type 2) 
intracellular signaling pathways. How the type 1/type 2 choice is made remains a central 
problem in the biology of apoptosis and has implications for apoptosis related diseases and 
therapy. We study the problem of type 1/type 2 choice in silico utilizing a kinetic Monte 
Carlo model of cell death signaling. Our results show that the type 1/type 2 choice is linked 
to deterministic versus stochastic cell death activation, elucidating a unique regulatory 
control of the apoptotic pathways. Consistent with previous findings, our results indicate 
that caspase 8 activation level is a key regulator of the choice between deterministic type 1 
and stochastic type 2 pathways, irrespective of cell types. Expression levels of signaling 
molecules downstream also regulate the type 1/type 2 choice. A simplified model of DISC 
clustering elucidates the mechanism of increased active caspase 8 generation and type 1 
activation in cancer cells having increased sensitivity to death receptor activation. We 
demonstrate that rapid deterministic activation of the type 1 pathway can selectively target 
such cancer cells, especially if XIAP is also inhibited; while inherent cell-to-cell variability 
would allow normal cells stay protected. 
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1. Introduction 

Apoptosis is the mode of cell death under a wide variety of cellular and physiological situations 
ranging from developmental regulations to tissue homeostasis. Dysregulation of apoptotic cell death 
has been implicated in a large number of diseases that includes degenerative disorders, cancer, and 
autoimmune diseases. It has long been known that apoptosis is coordinated through two distinct 
pathways: the type 1 (extrinsic) pathway that does not involve mitochondria, and the type 2 (intrinsic) 
pathway that requires mitochondrial activation. The extrinsic pathway is also called the death receptor 
mediated pathway as death receptor activation frequently induces apoptosis through this pathway. For a 
set of type 2 specific apoptotic stimuli (such as under DNA damaged conditions) only the type 2 
pathway is selectively activated. However, for a wide range of cellular and physiological situations, 
apoptosis is triggered by the activation of caspase 8 (initiator caspase) and thus could involve both 
pathways of apoptosis. Any type of death receptor (Fas, TNFR, DR4, DR5) mediated apoptosis falls 
into this category. A long-standing question in the biology of apoptosis is how the two pathways (type 1 
and type 2) are differentially activated [1-5]. Some of the initial studies indicated that caspase 8 
activation level coordinate the activation through those two pathways [1,6]. Strong caspase 8 activation 
was implicated in the type 1 choice whereas weak activation of caspase 8 was thought to choose the 
type 2 pathway for signal amplification. It was also thought that caspase 8 activation is cell type 
specific and cells were labeled as either type 1 or type 2 depending on their choice of the activation 
pathway. We have utilized Monte Carlo simulations to elucidate that for low caspase 8 concentration, 
the activation is dominated by the type 2 pathway with slow all-or-none type activation and large 
cell-to-cell variability, while for large caspase 8 concentration, the type 1 pathway is activated in a 
rapid deterministic manner [7]. These findings make the problem of type 1/type 2 choice even more 
intriguing as it becomes linked to deterministic/stochastic choice in apoptosis activation. The stochastic 
type 2 to deterministic type 1 transition was shown to be a robust feature of apoptosis signaling 
irrespective of cell types [8]. A recent experiment indicated that the expression of death receptors 
coordinate the type 1/type 2 choice and both pathways can be activated in any cell types by regulating 
the expression of death receptors [5]. Variation in the expression of death receptors presumably 
regulates the activation level of caspase 8 and thereby determines between type 1 and type 2 pathways. 
However, the mechanisms for death receptor mediated caspase 8 activation have not been clearly 
elucidated. It is also not clear how expression levels of various downstream signaling molecules in the 
two apoptotic pathways affect the type 1/type 2 choice. 

Previous studies in apoptotic cell death and apoptosis related diseases were dominated by population 
level measurements. More recent studies, however, indicate the presence of large cell-to-cell stochastic 
variability in the type 2 pathway of apoptosis [7,9-1 1]. It has been shown that cell-to-cell variability in 
type 2 apoptotic activation is characterized by slow progression but, eventual, all-or-none type 
activation for single cells. Inter-cellular variability in the expression levels of apoptotic proteins, such 
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as due to stochastic gene regulations, may introduce variability in the slow progression of apoptotic 
activation. This type of variability in apoptosis progression, coupled with rapid reaction events of 
cytochrome c release or caspase 9 activation, can lead to all-or-none type behavior. However, even 
when all the cellular parameters remain identical, inherent variability in apoptotic signaling reactions 
(due to small number of molecules or low-probability reactions) is capable of generating large 
cell-to-cell stochastic variability [7-10,12]. Both types of cell-to-cell variability can have important 
ramifications for diseases characterized by dysregulated apoptosis. In cancer, for example, cellular 
variation in anti-apoptotic protein levels (such as the Bcl2 concentration) as well as effective low 
number of molecules generated by over-expressed apoptotic inhibitors contribute to large cell-to-cell 
variability in apoptotic activation [10,13,14]. 

The programmed cell death mechanism of apoptosis utilizes a complex intracellular signaling 
network to carry out cellular demise in a controlled manner. In addition, phagocytic clearance of 
apoptotic bodies does not usually generate inflammatory conditions [15]. Targeting the apoptotic 
pathway is emerging as a new frontier in therapies of many of the apoptosis related diseases, such as 
cancer. However, large cell-to-cell variability, including inherent variability, can lead to fractional cell 
killing under chemotherapy and thus pose a challenge for developing cancer therapies that target the 
apoptotic pathway. Incomplete activation of the type 2 pathway may lead to generation of more 
resistant phenotypes. Therefore, it is important to search for strategies that would eliminate cell-to-cell 
variability in cancer cell apoptosis [13,14,16]. In addition, such strategies need to selectively target 
cancer cells. In this context, it is important to note that the expression levels of pro- and anti-apoptotic 
proteins vary significantly between normal and cancer cells [17-20]. In some cancer cells, which are 
equipped with increased levels of pro-apoptotic proteins Bid and Bax, it might be possible to reduce 
cell-to-cell variability by inhibiting Bcl2-like anti-apoptotic proteins and altering the Bcl2 to Bax ratio. 
Stochastic variability in Bax activation can even be abolished if a stochastic-to-deterministic transition 
is achieved [12]. However, in cancer cells for which the Bcl2 to Bax ratio has reached a very high level, 
alternative strategies such as switching the activation from type 2 to type 1 [7] may provide an option. 
Thus the question of type 1/type 2 choice is highly relevant for cancer cell apoptosis and cancer therapy. 

Mathematical and computational approaches, mainly based on ordinary differential equations and 
systems theoritic analysis, provided insight into the systems biology of apoptosis regulation [6,11,21-35]. 
Some of these studies have also considered death receptor mediated activation of the type 1 and type 2 
pathways. However, presence (or emergence) of small numbers of molecules (such as low cFLIP levels 
and small numbers of active caspase 9 molecules as a result of XIAP inhibition) makes consideration of 
stochastic effects essential in apoptosis modeling. Our Monte Carlo study of apoptosis signaling 
indicates that the type 2 pathway starts dominating when active caspase 8 concentration is low (~nM). 
Existence of low probability reaction events (such as the effective low rate of apoptosome formation) 
provides additional mechanisms for stochastic activation. 

In this work, we utilize Monte Carlo simulations to explore the systems level regulatory mechanisms 
of the type 1/type 2 choice. Expression levels of various proteins involved in the regulation of type 1 
and type 2 pathways were varied in isolation or in combination, and, activation of type 1 and type 2 
pathways were measured under those conditions. Results obtained from our in silico studies show that 
the type 1/type 2 choice is regulated at a systems level by coordinated expression levels of signaling 
molecules in apoptotic pathways. Concentration of active caspase 8 (initiator caspase) emerges as a key 
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regulator of the type 1/type 2 choice, consistent with previous studies [1,6,7]. Our results indicate a key 
role of the apoptotic inhibitor XIAP, as well as the XIAP to Smac ratio, in the type 1/type 2 choice and 
systems level regulation of apoptosis [3,4]. The formation rate of apoptosome is also shown to be 
important as its slow formation is a key rate limiting step in the type 2 pathway. 
In cancer cells, altered expression of various pro- and anti- apoptotic signaling proteins impact the 
type 1/type 2 choice. We demonstrate that increased sensitivity to death receptor activation in certain 
cancer cells can allow selective targeting of those cells (such as by death ligands) resulting in selective 
activation of caspase 8 in only those cells. XIAP inhibition in such death ligand treated cancer cells 
may result in a mixed type 1 -type 2 (or type 2) to type 1 transition in apoptotic activation and thus 
elimination of large cell-to-cell stochastic variability. 



2. Experimental Section 

2.1. The Signaling Model for Apoptotic Cell Death 

A detailed computational study is carried out utilizing kinetic Monte Carlo (MC) simulations of pre- 
and post-mitochondrial signaling events [7]. A simplified network model of apoptosis signaling is 
studied that is triggered by active capsase 8 (Figure 1) [4]. In some of the in silico experiments active 
caspase 8 was assumed to be present at initial time. To study apoptosis induction in cancer cells having 
heightened sensitivity to death receptor activation we incorporated a simplified model of caspase 8 
activation into our signaling model for type 1 and type 2 pathways. 

Figure 1. Schematic of the apoptotic death signaling network. Apoptosis is activated 
through two distinct pathways: type 1 (intrinsic) and type 2 (extrinsic). The type 1-type 2 
signaling loop is initiated by generation of active caspase 8 and ultimately converges on 
caspase 3/7 activation. 




Caspase 8 activation is known to be mediated by the clustering of adaptor proteins (such as 
FADD/TRADD) recruited to death receptor-ligand complexes. Procaspase 8 molecules are recruited to 
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the clustered adaptor proteins to generate the assembly of DISC (death-inducing-signaling-complex) 
and generate active caspase 8 molecules through autoprocessing [36]. In the current study, a simplified 
model of DISC (death inducing signaling complex) formation is considered where adaptor molecules 
can cluster (to lower thermodynamic free energy) when they are bound to death receptor-ligand 
complex; we call this state (receptor-ligand complex bound) of the adaptor molecule an active state. 
The parameter that captures the reduced energy of two neighboring adaptor molecules in active state is 
denoted by E D d (is taken to be -2 K B T unless specified otherwise). DISC formation is incorporated into 
the simulation by a hybrid simulation scheme between kinetic Monte Carlo model of intracellular 
signaling with an explicit free energy based model for the clustering of adaptor molecules [37,38]. 
Effective probability parameters P on and P 0 ff are introduced that capture an adaptor molecule's 
switching between an active and an inactive state (to capture the effect of death ligand induction such 
as FAS/TRAIL binding). Simulations are carried out for various values of the parameters P on and P 0 ff 
(those presumably vary depending on the cell type and/or the receptor type). 

Active caspase 8 initiates signaling through both type 1 and type 2 pathways. In the type 1 pathway, 
caspase 8 directly processes procaspase 3 to generate active caspase 3. In the type 2 pathway, caspase 8 
cleaves Bid to an active form (tBid) which translocates to mitochondria to bind with Bax. When two 
Bax molecules are bound to tBid (on the mitochondrial membrane) they could detach as an active Bax 
dimer. Apoptotic inhibitor Bcl2 molecules bind with tBid and Bax and thereby inhibit formation of 
active Bax dimers. It is also possible for Bid to directly activate Bax albeit with a low probability [12]. 
Cytochrome c is released into the cytosol in an all-or-none manner when the number of active Bax 
dimers reaches a pre-assigned threshold value [39,40]. Cytochrome c release leads to cytochorme 
c-Apaf binding and the subsequent formation of multi-molecular cyto c-Apaf-ATP complex 
apoptosome. Formation of the apoptosome complex is modeled in a simplified manner (by a cyto 
c-Apaf-Apaf-cyto c complex) where Apaf represents the Apaf-ATP complex. In Monte Carlo 
simulations, cytochrome c and Apaf molecules need to have spatial proximity before a binding reaction 
can occur between them, and such a diffusion limitation introduces additional probability into the 
cytochrome c-Apaf binding reaction. Low probability of apoptosome formation generates stochastic 
variability in apoptotic activation, but once formed, induces rapid activation of downstream caspases 9 
and 3 [7,41]. Procaspase 9 binds to the complex apoptosome and gets processed to its active caspase 9 
form, which in turn, cleaves procaspase 3 to its active form caspase 3. Apoptotic inhibitor XIAP binds 
to procaspase 9, caspase 9, and caspase 3, thus it could inhibit apoptotic activation in both type 1 and 
type 2 pathways. However, the release of mitochondrial Smac can antagonize the inhibition of XIAP. 
Smac is released simultaneously with cytochrome c in an all-or-none type manner. Appendix I contains 
all the reaction moves considered in our simulation of the apoptosis pathway. Activation of caspase 3 
(effector caspase) can be taken to be a downstream readout of apoptotic cell death signaling; MC 
simulations are carried out to measure the time-course of caspase 3 activation at a single cell level. 

Even though the present model can capture some of the key complexities of the apoptotic signaling 
pathway, it has several simplifying assumptions. For example, presence of multiple functionally similar 
proteins and their varied expression levels (depending on cell type) make the apoptotic pathway more 
complex than the present model can address. In our model, effects of functionally similar proteins are 
captured by a representative protein. Bcl-2 (B cell lymphoma protein 2) represents all the Bcl-2 family 
proteins (such as Bcl-2, Bcl-xL, Mcl-1) with similar anti-apoptotic properties. We do not explicitly 
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simulate BH3 only sensitizers (such as Bad or Bik), whose effect can be simulated by varying Bcl-2 
concentrations. However, coarse-graining the pathway by representative proteins allows us capture 
some of the essential biology of type 1/type 2 choice in apoptotic regulation. 

2.2. Monte Carlo Simulation of Cell Death Signaling 

Monte Carlo approach is well suited to simulate some of the complexities of signaling reactions such 
as the effect of spatial heterogeneity. Each run of our Monte Carlo simulation corresponds to a single cell 
activation, thus, Monte Carlo can capture cell-to-cell stochastic variability including inherent variability. 

We utilize a hybrid Monte Carlo simulation scheme that combines the following approaches: 
(1) a probabilistic rate constant based (implicit free energy) kinetic Monte Carlo simulation for various 
reaction moves such as diffusion, binding/unbinding and catalytic cleavage; (2) an explicit free-energy 
based model that captures clustering of death adaptor proteins (DISC clustering) utilizing 
energy-function based diffusion moves. Both of the above approaches have been utilized in previous 
works from this lab [38]. In our previous studies of apoptotic death signaling, a kinetic Monte Carlo 
based approach was taken which allowed straightforward estimation of probabilistic rate constants from 
the known diffusion and kinetic rate constants. Diffusion moves that are governed by attractive 
interactions with nearest neighbor molecules can be implemented in a simple manner through an 
explicit free energy based method (such as Lattice gas models in statistical mechanics [42]). Thus, an 
explicit free energy based model is utilized for simulating the diffusion moves of death adaptor 
proteins. At each Monte Carlo (MC) step molecules are randomly sampled 2N number of times, where 
N is the total number of molecules (either free or complexed) present in the system. Therefore, at each 
Monte Carlo (MC) step, one molecule is sampled (on average) twice to allow for one diffusion and one 
reaction move (only one move per molecule was allowed in our previous works). All the reaction 
moves considered in our simulations are provided in Appendix I (corresponding kinetic rate constants 
are also included). Diffusion moves are carried out to one of the randomly chosen neighboring sites 
(4 for membrane bound molecules and 6 for cytosolic molecules) provided the condition of mutual 
physical exclusion is satisfied. In the kinetic Monte Carlo part of the simulation, once a 
diffusion/reaction move is randomly sampled it is accepted only if a randomly generated number in 
[0,1] is less than the pre-defined probability constant for that particular move (otherwise the move is 
rejected). For reversible reactions, such as the binding/unbinding reactions A + B ^± C, the detailed 

balance condition is satisfied by the choice of P on and P 0 ff (through the ratio P on /Poff), the probability 
constants for the binding and unbinding reactions, respectively (see [38] for details). In the explicit free 
energy based part of the simulation, diffusion moves are accepted based on Metropolis criterion 
(Paccpet = min[l,exp(-AE/K B T)], where AE is the free energy difference between the new and the 
previous configuration [42]). 

2.3. Estimation of Parameter Values (as Used in Our Monte Carlo Simulation) 

A simulation volume of 1.2 x 1.2 x 1,2 um 3 (corresponding to a 60 x 60 x 60 lattice with lattice 
spacing Ax -20 nm) is chosen in such a manner that the number of molecules (for each molecular 
species) is equal to the nanomolar concentration. Cytochrome c/Smac is initially contained in a 
mitochondrial volume of 0.36 x 0.36 x 0.36 um (18 x 18 x 18 lattice). Utilization of a small system 
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size significantly reduces the computational cost of the simulation. Each MC step (AT) is chosen to be 
2 x 10~ 4 s based on known mobility of cytosolic molecules. This allows us to take the probability of 
diffusion (Pdiff) for the fastest diffusing species (cytosolic molecules) to a value 0.1 (an approximate 

2 2 

diffusion constant can be estimated D ~ 0.1 x (Ax) /(AT) ~ 0.2 \im /s). Consistent with lower diffusion 
constants known for membrane proteins, probability of diffusion for molecules on the plasma 
membrane (or the mitochondrial membrane) is taken an order of magnitude lower than that is used for 
cytosolic molecules. It is reasonable to expect that the multi-molecular complex apoptosome will have 
significantly reduced mobility and its Pdiff is assumed to be zero. Kinetic reaction rates (such as k on /k 0 ff) 
and molecular concentrations (Tables Al and A2 in Appendix I) are obtained from values reported in 
the literature and utilized in our previous work [7,9,10,12,13] (unless specified otherwise). Probabilistic 
reaction rate constants used in this study are obtained from kinetic rate constants using a previously 
described parameter-mapping scheme [7,38]. P 0 /f (or P cat ) simulation parameters are obtained by 
multiplying k 0 /f (or k cat ) values by 10" s (0.5 MC time-step). Probabilistic parameters for association 
reactions are determined using the following relation: P on = 10 x k on nM -s [7,38]. Values for 
probabilistic simulation parameters are provided in Table A3 (Appendix I). 

A typical simulation is run until the concentration of unbound active caspase 3, the effector caspase, 
reaches its half-maximal value ~50 nM (simulation stops if it is not completed before a predefined 
number of MC steps). Certain threshold amount of caspase 3 activation might be sufficient for 
downstream PARP cleavage and that threshold should also regulate the type 1/type 2 choice. In this 
work, data normalized to its half-maximal value is reported for caspase 3 activation. Type 1/type 2 
classification is also based on half-maximal activation of caspase 3. Controlled Monte Carlo 
experiments are carried out for specific parameter values (such as molecular concentrations). Each run 
of the simulation corresponds to activation at a single cell level. 

3. Results and Discussion 

3.1. Increased Active Caspase 8 Concentration Switches the Activation from Stochastic Type 2 to 
Deterministic Type 1 

We have previously shown that two apoptotic pathways are activated in a distinct manner [7,8]. 
Apoptotic cell death can be activated in rapid deterministic manner through the type 1 pathway. 
In contrast, apoptotic activation through the type 2 pathway is usually slow with inherent cell-to-cell 
variability. How these two distinct pathways are activated remains an unresolved issue in the biology of 
apoptosis. Previous studies indicated the role of active caspase 8 as a key determinant of the type 1/type 2 
choice [1,2,6,7]. In Figure 1, the signaling network for apoptotic cell death is shown where the type 1 
and type 2 pathways form a loop network that can be triggered by active caspase 8. Clearly, caspase 8 
is expected to be a critical regulator of the type 1/type 2 choice and such a task is achieved by 
differential binding to its immediate binding partners procaspase 3 (type 1 pathway) and Bid (type 2 
pathway). Capsase 8 binds to Bid with moderately high affinity (k A = k on /k 0 ff ~ 10 9 M -1 ) compared to 
its relatively weak affinity for procapsase 3 (k A = k on /k 0 ff ~ 1.67 x 10 5 M _1 ). Thus for low 
concentrations of active caspase 8 (weak apoptotic stimuli) the type 2 pathway is preferentially 
activated (Figure 2a). Large caspase 8 concentration (strong apoptotic stimuli), on the other hand, is 
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sufficient to directly activate procaspase 3 in the type 1 pathway (Figure 2b, c). Mixed type 1-type 2 
behavior is also observed (where procaspase 3 is cleaved by both caspase 8 and caspase 9), especially 
for intermediate levels (~5 nM) of caspase 8 activation (Figure 2b). A large contribution to cell-to-cell 
variability in caspase 3 activation (Figure 2a) seems to originate from stochastic variability in 
apoptosome formation (post-mitochondrial signaling module) but remains sensitive to its rate of 
formation [9]. The type 2 pathway is designed to amplify an initially weak signal (to a strong 
all-or-none type activation) while large cell-to-cell variability generates heterogeneity in cellular 
response presumably as an adaptive strategy to a weak stimulus [7,8]. 

Figure 2. Caspase 3 activation for three different caspase 8 concentrations: 1 nM (a), 5 nM 
(b), 10 nM (c). Data is shown for 10 representative single cell runs; each color corresponds 
to apoptosis activation for a single cell (Monte Carlo run). Type 2 activation is 
characterized by all-or-none type behavior with large cell-to-cell variability. Increasing 
concentration of caspase 8 results in increased type 1 activation and suppression of 
cell-to-cell variability. 




time (MC steps) 

In the limit of weak capsase 8 activation (~1 nM), a set of ODEs describing caspase 8 reactions to 
Bid and procaspase 3 can be utilized to approximately estimate that at early times caspase 8-Bid 
complex is significantly more abundant than caspase 8- procaspase 3 complex (Appendix II). It can 
also be shown that (at early-times) unless caspase 8 concentration is very high (-100 nm or higher), 
caspase 8 starts selectively activating the type 2 pathway by cleaving Bid to tBid. Even though the 
initial type 1/type 2 preference is determined by the amount of active caspase 8, the final type 1/type 2 
choice should depend on how fast activation signal through the two pathways can lead to completion of 
caspase 3 activation. 
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3.2. Systems Level Regulatory Mechanisms of Type 1 /Type 2 Choice 

Our previous studies elucidated mechanisms for slow activation and generation of cell-to-cell 
variability through the type 2 pathway: (1) effective small number of reactant molecules generated due 
to the action of inhibitory molecules such as Bcl2 and XIAP [10,13] and (2) low probability reaction 
events such as the apoptosome formation [7]. These parameters are varied here to study their impact on 
the type 1/type 2 choice. 

3.2.1. Effect of Small (~3-Fold) Bcl2 Variation on Cyto C Release and Type 1/Type 2 Choice 

Small variations in Bcl2 levels can capture tissue specific variations as well as variations during 
various developmental stages [43]. Bcl2 over-expression (3-fold) resulted in inhibition of Bax 
activation and slightly delayed release of mitochondrial cytochrome c. Subsequent formation of 
XIAP-Smac and apoptosome complex were also delayed. For caspase 8 = 1 nM, the average time to 
apoptosome formation was higher when Bcl2 was 3 -fold overexpressed compared with the same for 
regular Bcl2 expressions (2.74 x 10 7 vs. 2.39 x 10 7 MC steps). Slower apoptosome formation resulted 
in a slightly delayed type 2 activation and a small increase in type 1 activation. When active caspase 8 
concentration was increased to 10 nM, type 1 pathway dominated the activation and Bcl2 inhibition no 
longer had a significant effect. Bcl2 is a known oncogene and large over-expression (> 10-fold) of Bcl2 
is considered later in the context of cancer cells. 

3.2.2. XIAP to Smac Ratio Is a Key Regulator of Type 1/Type 2 Choice: Inhibition of the Type 1 
Activation by XIAP 

XIAP is known to be a key regulator of the apoptotic activation as it binds and inhibits both 
caspase 9 (including procaspase 9) and caspase 3. Previous experimental studies indicated that under 
XIAP inhibition caspase 8 processing of procaspase 3 cannot proceed to its fully cleaved pl7/pl2 form, 
instead it remains in a partially processed p20/pl2 form [3]. It was observed that XIAP could inhibit 
both pl7/pl2 and p20/pl2 forms of active caspase 3. It was also shown that mitochondrial release of 
Smac, upon threshold Bax activation, can inhibit XIAP and thus modulates the type 1/type 2 choice in 
apoptosis activation [3]. We carried out simulations for various XIAP/Smac ratios [3,44]. Initially 
XIAP concentration was kept constant at 30 nM while the Smac level was varied: 50 nM (XIAP to 
Smac ratio 0.6) and 10 nm (XIAP to Smac ratio 3). When caspase 8 = 1 nM, activation is dominated by 
the type 2 pathway as observed in type 2 cells. Decreasing the Smac concentration from 50 to 10 nM 
resulted in increased XIAP inhibition of caspases and suppression of the initial type 1 activation 
(Figure 3a,b; Table 1). Similar behavior was observed when XIAP concentration was increased to 90 nM 
(Smac = 50 nM) (Figure 3c; Table 1). In this case, the type 2 activation was delayed (Figure 3c) until 
the XIAP-caspase 3 complexes reached sufficiently high concentration so that procaspase 9 molecules 
could be relieved from XIAP inhibition. For XIAP = 60 nM and Smac = 50 nM, removing Bid from the 
system resulted in weak type 1 activation (Figure 4a). However, when XIAP was inhibited (in Bid 
deficient cells), slow type 1 activation was observed (Figure 4b, c) consistent with earlier experimental 
studies of Fas ligand induced apoptosis in Bid -/ ~XIAP~ /- mice hepatocytes (type 2 cells) [4]. When the 
initial XIAP concentration is set to zero in presence of Bid, more rapid apoptotic activation was 
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observed and for a few cells half-maximal activation was nearly completed through the type 1 pathway 
(Figure 5). Type 1 activation became increasingly dominant for increased active caspase 8 (-2-3 nM) 
(Figure 5). Even in presence of significant amount of XIAP, when active caspase 8 concentration was 
increased to ~5 nM (or higher) significant type 1 activation was observed along with the type 2 
signaling (Figure 3d,e; Figure 6). Mitochondrial Smac release inhibited XIAP and assisted type 1 
activation to proceed. Such a mode of apoptotic activation can be labeled as type 2 assisted type 1 
activation. Eventually, caspase 9 was activated and subsequent activation of caspase 3 progressed 
through both type 1 and type 2 pathways (mixed type 1-type 2 behavior) (Figure 6). When XIAP was 
highly expressed (-90 nM), type 1 activation was observed to proceed only after XIAP was inhibited 
by mitochondrial Smac release (Figure 3f) [3]. The effect of type 1 assisted type 2 activation was 
diminished when Smac release was suppressed in a significant manner. In addition, reduced Smac 
concentration resulted in slower rate of caspase 9 activation (as more XIAP becomes available to 
inhibit procaspase 9) and slower apoptosis. Therefore, apoptotic inhibitor protein XIAP (also the XIAP 
to Smac ratio) emerges as a key determinant of the type 1/type 2 choice in apoptotic activation. The 
impact of XIAP inhibition in type 1/type 2 choice can be explained by noting that XIAP inhibits several 
apoptotic molecules in both type 1 and type 2 pathways; in addition, it inhibits the effector caspases 
(such as caspase 3) that close the type 1-type 2 loop network. 

Figure 3. Single cell caspase 3 activation for different values of XIAP to Smac ratio and 
active caspase 8 levels. Two distinct active caspase 8 concentrations are studied: 1 nM (left 
panel: a,b,c) and 10 nM (right panel d,e,f). The following XIAP and Smac levels are 
simulated: XIAP = 30 nM, Smac = 50 nM (a and d), XIAP = 30 nM, Smac = 10 nM (b and e), 
XIAP = 90 nM, Smac = 50 nM (c and f). Data is shown for 5 representative single cells in 
each of the above cases; each color corresponds to apoptosis activation for a single cell 
(Monte Carlo run). 




5e7 

time (MC steps) 



5e7 

time (MC steps) 



1e8 
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Table 1. Caspase 3 activation for different XIAP/Smac ratios and at different time -points. 
Active caspase 8=1 nM. Data normalized to half-maximal caspase 3 (= 50 nM) activation 
is shown. Data is averaged over 64 single cells (MC runs) for XI AP = 30 nM and over 10 
single cell runs for XIAP = 90 nM. 



XIAP/Smac 


Average Caspase 3 activation (SD) 




combinations 


T = 2.5 x 10 7 (MC steps) 


T = 5 * 10 7 (MC steps) T = 


1 x 10 8 (MC steps) 


XIAP = 30 nM 


0.65 (0.38) 


0.92 (0.20) 


1.0 (0.0) 


Smac = 50 nM 








(ratio 0.6) 








XIAP = 30 nM 


0.47 (0.38) 


0.81 (0.32) 


0.98 (0.08) 


Smac = 10 nM 








(ratio 3) 








XIAP = 90 nM 


0.01 (0.01) 


0.1 (0.08) 


1.0 (0.0) 


Smac = 50 nM 








(ratio 1.8) 









Figure 4. Single cell caspase 3 activation in Bid deficient cells for the following values of 
XIAP inhibition: XIAP = 60 nM, (a), XIAP = 30 nM, (b), XIAP = 0 nM (c). Smac 
concentration is kept constant at 50 nM. Data is shown for 10 representative single cells for 
each XIAP concentrations; each color corresponds to apoptosis activation for a single cell 
(Monte Carlo run). Decreasing XIAP concentration results in increased type 1 activation. 




time (MC steps) 
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Figure 5. Single cell caspase 3 activation in cells where XIAP inhibition is removed. Data 
is shown for 10 representative single cells for each of the following caspase 8 
concentrations = 1 nM (a), 2 nM (b) and 3 nM (c); each color corresponds to apoptosis 
activation for a single cell (Monte Carlo run). 




time (MC steps) 

3.2.3. Increased Rate of Apoptosome Formation Favors the Type 2 Pathway 

To study the effect of slow apoptosome formation on type 1/type 2 choice, we varied the Apaf 
concentration (-100 nM in reference [6] and ~20 nM in reference [26]; [45,46]) and the kinetic reaction 
rate k on for the cytochrome c — Apaf association reaction (k on may vary under pH change [47]). A key 
role of apoptosome formation on type 2 signaling has also been noted in an earlier study [25]. For low 
concentrations of caspase 8 (=1 nM) and when Apaf = 100 nM, activation is dominated by the type 2 
pathway (Figure 2a). When Apaf concentration was reduced to 20 nM, slow type 1 activation 
frequently replaced type 2 activation (Figure 7a). For low Apaf concentrations, low-probability 
apoptosome formation becomes very slow, allowing type 1 activation to complete before caspase 9 is 
activated (Figure 7). When the low probability rate constant of cytochrome c and Apaf binding was 
assumed 2-fold higher, formation of apoptosomes was faster, leading to increased type 2 activation but 
diminished cell-to-cell variability. In contrast, when the cyto c-Apaf binding rate was taken 2-fold 
lower, increased type 1 signaling was observed due to slower rate of formation of apoptosomes. The 
average times to form the first apoptosome are provided in Table 2. For a large (~ 10-fold or higher) 
increase in the cytochrome c-apaf association constant, cell-to-cell variability generated through the 
post-mitochondrial signaling module is significantly decreased. 
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Figure 6. Fraction of apoptotic cells (measured by half-maximal caspase 3 activation) in 
WT and Bid _/ " cells at 2.5 x 10 7 (a), 3.75 x 10 7 (b) and 5 x 10 7 (c) MC steps. Data is 
obtained from 64 single cell experiments (Monte Carlo runs) with active caspase 8-5 nM. 
Such moderate level of caspase 8 activation might be possible to achieve in type 1 cells (at 
least for early times). When the type 2 pathway is blocked (by removing Bid), strong 
suppression of caspase 3 activation was observed for early times but late time apoptosis was 
not inhibited. (Results from our in silico studies can be compared with previous 
experimental data in type 1 thymocyte cells [4]: ICAD (inhibitor of caspase-activated 
DNase) cleavage and fluorogenic measurement of effector caspase activity (DEVDase) 
indicate suppression of apoptotic activation at early time (~1 h) in Bid deficient thymocytes. 
The early-time suppression of apoptosis could be due to lack of capsase 9 activation; it is 
also possible that reduced apoptosis in Bid deficient cells results from lack of type 2 
mediated inhibition of XIAP (by Smac/DIABLO). In our in silico experiments, both XIAP 
inhibition and caspase 9 activation contribute to early-time activation of effector caspases. 
Therefore, this data indicates a role of the type 2 pathway in apoptotic activation in type 1 
cells such as thymocytes). 




time = 2.5e7 MC steps 



time = 3.75e7 MC steps 



time = 5e7 MC steps 



Figure 7. Single cell caspase 3 activation with concentration of Apaf reduced to 20 nM (a). 
Data is shown for 8 representative single cells; each color corresponds to apoptosis 
activation for a single cell (Monte Carlo run). Caspase 9 activation (shown in b) is slower 
(compared with the same for Apaf = 100 nM) due to decreased apoptosome formation and 
allows type 1 activation to progress longer. 
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Table 2. Average time to first apoptosome formation as the rate constant for the 
cytochrome c-Apaf binding is varied. Data is averaged over 56 single cells (MC runs). 



Kinetic rate constant for Cyto 


Time to (first) apoptosome 


c-Apaf binding 


formation (MC steps) 


2-fold high 


1.03 x 10 7 


normal 


2.47 x 10 7 


2-fold low 


3.90 x 10 7 * 



* considering cells in which apoptosome has formed before caspase 3 activation reached its half maximal 
value (the actual average is higher than reported). 

3.3. Robustness of Type 1/Type 2 Choice in Apoptosis Regulation 

Previously we utilized a minimal signaling network (Figure 8) to demonstrate that stochastic 
type 2 to deterministic type 1 transition in apoptosis activation can be achieved irrespective of cells 
types (provided a few simple conditions are satisfied) [8]. The same minimal network is studied here to 
show that some of the systems levels regulatory mechanisms of the type 1/type 2 choice, as observed in 
Monte Carlo simulations of the present work, are robust features of apoptosis signaling (supplemental 
material). Specifically, we varied the following variables in the minimal network: (i) the concentration 
of the signaling molecule that opens the type 1/type 2 loop and (a measure of strength of the apoptotic 
stimulus) (ii) the rate constant for the slow step in the type 2 pathway. The signaling molecule that 
opens the type 1/type 2 loop in the minimal network captures the effect of active caspase 8 in the full 
apoptotic network. Increasing (decreasing) its concentration resulted in increased (decreased) activation 
of the type 1 pathway (Figure 9); result was shown to be robust over 2-fold variations in concentrations 
of other signaling molecules. Increasing (decreasing) the rate constant for the slow step in the type 2 
pathway, which could capture variation in the formation rate of apoptosome, enhanced (diminished) the 
type 2 activation (Figure 10). These results are consistent with findings of the in silico experiments of 
the apoptotic death pathway and elucidate some of the robust regulatory mechanisms of type 1/type 2 
choice in apoptosis signaling. The precise concentrations of different molecules (in the full apoptotic 
network) at which the type 1/type 2 transition occurs should depend on cell type specific features and 
can be better predicted as more accurate quantitative information (such as regarding molecular 
concentrations and rate constants) becomes available. 




Cells 2013, 2 



375 



Figure 8. A previously developed model of minimal network [8] is utilized to elucidate 
some of the robust systems level regulatory mechanisms in apoptosis signaling. This 
signaling network can be thought of as a highly coarse-grained representation of the full 
apoptotic signaling network. The following rate constants were used: 2 x 10~ 4 (for the type 1 
pathway); 1.0, 10~ 7 and 10" 1 (for the type 2 pathway). X 0 measures the strength of an 
apoptotic stimulus (varied in Figure 9). Initial numbers for all other the molecules (Xi, X2, 
X3) were assumed to be 100. 




Figure 9. The effect of concentration (X 0 ) variation of the molecule that opens the type 1-type 
2 loop network. Two different values ofX 0 were studied: 1 (left panel) and 20 (right panel). 
Results are shown to be robust for a 2-fold variation in concentration fluctuations for other 
molecules: normal concentration (a and d), 2-fold underexpression (b and e), 2-fold 
overexpression (c and f). The minimal network was analyzed using the stochastic differential 
equations developed in [2]. The SDE that captures low probability reaction of X 2 —> X 2 * 
(slow step) was solved using a Poisson Runge-Kutta scheme [12]. In this method, the number 
of times a specific reaction channel fires in a given time t is a Poisson random variable with 
mean At and variance At, where A is the propensity function. For the slow reaction under 
consideration A = k(dt)Xj*X2. The rate constant for the slow step is k and dt is the time-step 
used in the numerical solution. All other reaction SDEs were solved using the standard 
Euler-Maruyama numerical scheme for solving stochastic differential equations [8,12]. 
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Figure 10. The effect of variation in the rate constant for the slow step (JG —> X2*) is 
studied for the following values: 10~ 6 (a), 10~ 7 (b), and 10~ 8 (c). The minimal network was 
analyzed using the stochastic differential equations developed in [2]. The method utilized to 
solve the SDEs is the same as in Figure 9. X 0 (apoptotic stimulus strength) was taken as 1. 
(A similar study using Gillespie's stochastic simulation algorithm (SSA) was carried out [8]). 
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3.4. Regulation of Type 1 /Type 2 Choice in Cancer Cells Having Over-Expressed Anti-Apoptotic 
Proteins Bcl2 and XIAP 

A hallmark of cancer is genomic instability and aberrant expression of oncogenes and proteins. 
Expression levels of signaling proteins in the apoptotic death pathway are particularly altered, including 
the level of well known oncogenic proteins such as Bcl2, resulting in dysregulated death signaling and 
inhibition of apoptosis. Previous studies elucidated the mechanism of Bcl-2 and XIAP inhibition of 
type 2 apoptosis in cancer cells [13]. It was shown that increased expression of Bcl2 and/or XIAP 
enhances cell-to-cell stochastic variability and time-to-death in type 2 activation, imparting apoptosis 
resistance to cancer cells. In this work, we focus on the regulation of type 1/type 2 choice in cancer 
cells that are equipped with increased levels of Bcl2 and XIAP. We vary the concentrations of Bcl2 and 
XIAP and study type 1/type 2 control of apoptosis. A 2-fold over-expression for Bid and Bax molecule 
is also assumed [12]. 

Bcl-2 over-expression resulted in slower cytochrome c/Smac release with increased cell-to-cell 
stochastic variability, consistent with observations of previous studies [10,12]. When Caspase 8=1 nM 
and XIAP = 30 nM, 15-fold Bcl2 over-expression allowed initial slow type 1 activation to proceed 
longer but eventually switched to all-or-none type 2 activation (Figure 11a). Cell-to-cell variability in 
caspase 3 activation (Figure 11a) results from stochastic variability in (i) Cytochrome c/Smac release 
(pre-mitochondrial signaling module) and (ii) apoptosome formation (post-mitochondrial signaling 
module) (Table 3). When XIAP is 3-fold over-expressed, strong suppression of the initial type 1 
activation was observed and caspase 9 activation precedes caspase 3 activation (Figure lib). For high 
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XIAP over-expression (~5-fold of normal 30 nM), apoptotic activation is strongly suppressed (Figure 1 lc) 
providing a mechanism for resistant cancer cells. Strong suppression of apoptotic activation under high 
XIAP levels has also been observed earlier [27]. Certain cancer cells are known to be chemosensitive to 
death ligand induction presumably by having increased sensitivity to death receptor activation [19]. In 
such cancer cells it is expected that increased death receptor activation would result in increased 
amount of active caspase 8 generation. When caspase 8 = 10 nM but XIAP level remains normal 
(30 nM), rapid deterministic activation of the type 1 pathway is observed (irrespective of Bcl2 
over-expression levels) (15-fold Bcl2 in Figure lid). Therefore, achieving such a caspase 8 activation 
level (-10 nM) combined with lower XIAP level (-30 nM) can be an optimal strategy for targeting 
certain cancer cells. Increased XIAP expression, however, inhibits the type 1 pathway of apoptosis. 
When XIAP is 3-fold over-expressed along with 15-fold over-expression of Bcl2, type 1 activation is 
suppressed until Smac is released from mitochondria. Smac binding to XIAP allows caspase 3 
activation to proceed and subsequent activation of caspase 9 induces additional caspase 3 activation 
through the type 2 pathway (Figure lie). Hence, a type 2 assisted type 1 mode of activation followed 
by a mixed typel-type2 activation was observed when XIAP was moderately over-expressed (~3-fold). 
When Bcl2 over-expression was taken to be 50-fold (keeping XIAP ~3-fold), Cyto c/Smac release 
became slower with increased cell-to-cell stochastic variability. As a result, caspase 3 activation 
showed large cell-to-cell stochastic variability with all-or-none type behavior (Figure 12). Increasing 
the XIAP over-expression level to 5-fold led to strong XIAP mediated inhibition of caspase 3 activation 
and weak type 1 activation (Figure 12c,f). A detailed analysis of caspase 3 activation for various 
combinations of Bcl2 and XIAP overexpression is provided in Figure 13. Additional interactions such as 
caspase 6 mediated activation of caspase 8 (positive feedback) [3] and/or XIAP degradation of active 
caspase 3 can further impact the type 1/type 2 choice in a cell type dependent manner. 

Figure 11. Time course of caspase 3 activation in cancer cells having over-expressed Bcl2 
and XIAP. Bcl2 over-expression is taken as ~T5-fold while XIAP expression is varied: 
normal (a and d), 3-fold (b and e), 5-fold (c and f). Two different caspase 8 concentrations 
are used: 1 nM (left panel: a, b, c) and 10 nM (right panel: d, e, f). Data is shown for five 
representative single cells in each of the above cases; each color corresponds to apoptosis 
activation for a single cell (Monte Carlo run). 
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Table 3. Cell-to-cell stochastic variability in time-to-death (capsase 3 activation) results 
from cell-to-cell variability in (i) cytochrome c release (pre-mitochondrial signaling module) 
and (ii) apoptosome formation (post-mitochondrial signaling module). Time-to-Cytochrome 
c/Smac release is estimated from the initiation of XIAP-Smac complex formation. 
Time-to-apoptosome formation indicates the time-scale of first apoptotosome formation. 
Time-to-death is measured by the time by which active capsase 3 concentration reaches 50 nM 
(half-maximal level). Data is shown for three single cells (MC runs). Time is measured in 
MC steps. 





Cell-to-cell variability in caspase 3 activation 






Time-to-Cyto c /Smac release Time-to-apoptosome formation 


Time-to-death 


Cell 1 


6.2 x 10 6 8.6 x 10 6 


1.6 x 10 7 


Cell 2 


7.6 x 10 6 9.1 x 10 7 


9.5 x 10 7 


Cell 3 


7.0 x 10 6 4.2 x io 7 


4.7 x 10 7 



Figure 12. Time course of caspase 3 activation in cancer cells having over-expressed Bcl2 
and XIAP. Bcl2 over-expression is taken as ~50-fold while XIAP expression is varied: 
normal (a and d), 3-fold (b and e), 5-fold (c and f). Such high Bcl2 expression can result 
from combined effect of various Bcl2 like anti-apoptotic proteins (also possibly found in 
cancer stem cells). Two different caspase 8 concentrations are used: 1 nM (left panel: a, b, c) 
and 10 nM (right panel: d, e, f). Data is shown for five representative single cells in each of the 
above cases; each color corresponds to apoptosis activation for a single cell (Monte Carlo run). 
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Figure 13. Activated caspase 3 (fraction of half-maximal activation) for various 
combinations of Bcl2 (15- and 50-fold) and XIAP overexpression (normal, 3-fold and 
5-fold). Two different values of initial caspase 8 concentration are studied: 1 nM (top) and 
10 nM (bottom). Data obtained from 60-64 single cells (MC runs) are shown for two 
different time instants: 5 x 10 7 (a) and 1 x 10 8 (b) MC steps. 
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3.5. Apoptosis in Cancer Cells Having Heightened Sensitivity to Death Receptor Activation: A 
Mechanism for Selective Targeting of Cancer Cells 

Targeting the apoptotic pathway can be an effective way to eliminate cancer cells in a controlled 
manner. However, one needs to find a strategy that will selectively target cancer cells leaving normal 
cells protected. It is known that in certain cancer cells significant expression of death receptors (such as 
Fas/TNFR/TRAIL (DR4/DR5)) provides an opportunity to selectively target cancer cells [19,48,49]. 
In addition, normal cells might have protection from a larger expression of decoy receptors (compared 
to that in cancer cells) [50,51] resulting in the effective reduced expression of death receptors. The 
specificity for cancer cell apoptosis can be induced (or enhanced) by additional mechanisms such as 
gene therapeutic approaches [52]. The specific choice of death ligand (or agonist antibody) depends on 
a variety of factors such as whether the targeted death receptor generates type 1 or type 2 activation in 
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healthy cells [19]. Even though some aspects of activation and downstream effects may vary depending 
on the specific death receptor under consideration, the mechanism of death receptor mediated caspase 8 
activation (from procaspase 8) seems to be similar in many of those cases and involves (i) recruitment 
of adaptor molecules (to ligand bound death receptors) containing both death domain and death effector 
domain (such as FADD) and (ii) formation and clustering of DISC [36]. To capture selective killing of 
cancer cells, we study apoptosis mediated by death receptor induced caspase 8 activation. The 
mechanism of death receptor clustering (such as by death ligand induction) is modeled in a simplified 
manner by DISC clustering (see Experimental Section). To capture the effect of death ligand induction 
(such as by TRAIL), we carried out simulations for the various P on and P 0 ff values. Even though some 
clustering was started to be seen for P off values ~10~ 2 (P on =1), significant clustering was observed 
when P 0 ff = 10~ 3 (such as in cancer cells). DISC clustering is shown in the supplemental material 
(Figure 14). Corresponding caspase 8 activation is shown in Figure 15. Lower values for the DISC free 
energy parameter Edd resulted in increased caspase 8 activation (Figure 15). The average time-to-death 
(Td) decreased with increased DISC clustering and caspase 8 activation (Td = 4.46 x 10 7 MC steps for 
Edd = 0; T d = 2.95 x 10 7 MC steps for E dd = -K B T; T d = 2.93 x 10 7 MC steps for E dd = -2K B T). Slower 
generation of active caspase 8 also provides a mechanism for increased cell-to-cell stochastic 
variability. For normal cells, reduced expression of death receptors and the inhibitory effect of decoy 
receptors is captured by lower P on values of 10~ 4 (16% cell death at t = 10 8 MC steps) and 10~ 5 
(1.6% cell death at t = 10 8 MC steps) and P off = 1. 

Figure 14. DISC (death inducing signaling complex) clustering for three different values of 
the DISC clustering energy parameter: E d d = 0, -K B T and -2K B T (a, b, c). Adapter 
molecules (that bind death-receptor ligand complex to generate DISC) are shown on the cell 
surface of 60 x 60 lattice sites (1.2 jam x 1.2 jam). Pon = 1 and Poff = 10~ 3 are used in the 
simulations. Representative single cell data are shown after 10 7 MC steps. 
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Figure 15. Activated caspase 8 (fraction of maximal) for three different values of DISC 
clustering energy: (a) Edd = 0, (b) Edd = _ K B T, (c) Edd = -2K B T. Pon = 1 and Poff = 10" 
are assumed. We simulated cancer cells having over-expressed Bid (2-fold), Bax (2-fold), 
Bcl2 (15-fold) and XIAP (3 -fold). Data analyzed from 64 single cells (MC runs) are shown 
for three different time instants: 5 x 10 6 , 10 7 and 3 x 10 7 MC steps. Clustering of DISC 
(when Edd = ~K B T or -2K B T) resulted in increased active caspase 8 generation and 
subsequent activation of effector caspases. 
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The mechanism of death receptor clustering (such as by death ligand induction) and subsequent 
caspase 8 activation is complex. Formation of closed oligomers (such as of Fas-FADD [53]) within the 
DISC and/or altered partitioning of membrane proximal molecules in lipid rafts (upon ligand induction) 
may lead to a stronger effect of DISC formation on caspase 8 activation [54,55]. A detailed model that 
will explicitly simulate death receptor-ligand binding and lipid mediated interactions is currently under 
development. The strength of an external apoptotic stimulus, such as concentrations of a death ligand, 
is expected to affect death receptor activation and the type 1/type 2 choice, but the mechanism remains 
to be explored. Both apoptotic activation as well as protection (from apoptosis) due to over-expressed 
Bc1-xl has been shown to correlate with death ligand concentrations [3,5]. Studies in neural cells 
indicated activation of both pathways when treated with A[3-oligomers (type 2 activation) and 
A^-aggregates (type 1 activation) [56]. In addition to its key role in apoptosis, mechanisms of death 
receptor activation and caspase 8 generation have implications for other modes of cell death such as in 
necroptosis [57,58]. 

Cancer drugs that target the apoptotic pathway by perturbing the membrane lipid environment have 
attracted recent attention [59,60]. The mechanisms of death receptor activation may also have 
implications for the cancer stem cell model. Recent studies in severely immunocompromised 
NOD/SCID IL-2rg _/ ~ mice found significantly increased fraction of tumorigenic cells (melanoma stem 
cells) when compared with the same in NOD/SCID mice [61]. Cellular and molecular basis of such 
increased fraction of cancer stem cells in NOD/SCID IL-2rg _/ ~ mice is yet to be elucidated. However, a 
role of NK cell mediated tumor cell killing can be postulated as NOD/SCID IL-2rg~ /_ mice are known 
to be depleted of NK cells (due to the null Interleuken-2 receptor gamma chain mutation). NK cells are 
equipped with death ligands that are known to activate the death receptors DR4 and DR5 on target 
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cells [62]. Therefore, NK cell mediated apoptotic activation may provide a mechanism for increased 
killing of tumor cells transplanted in NOD/SCID mice. 

3.6. Combined Death Ligand Induction and XIAP Inhibition can Be an Optimal Strategy to Kill Cancer 
Cells: Maximizing Specificity and Minimizing Cell-to-Cell Stochastic Variability 

Highly over-expressed anti-apoptotic proteins allow cancer cells to resist apoptotic cell death. 
Increased inherent cell-to-cell variability, in apoptotic activation of cancer cells equipped with 
over-expressed anti-apoptotic proteins, can be a mechanism for the resistant phenotype of cancer cells. 
Fractional killing of cancer cells (under chemotherapy) need to be addressed as survived cells may 
acquire more resistant phenotypes. In a recent work we have shown that over-expression of certain 
pro-apoptotic proteins, along with over-expressed apoptotic inhibitors, can allow a stochastic-to- 
deterministic transition in the type 2 pathway (by simply inhibiting Bcl2 like proteins in the type 2 
pathway) [12]. While such a strategy might be effective for many types of cancer cells, there are cancer 
cells that will remain resistant due to very high Bcl2 to Bax ratio. For cancer cells that are resistant to 
Bcl2 inhibitor chemotherapy, switching the activation from type 2 to type 1 can be a potential option 
for inducing rapid deterministic apoptosis. Such a strategy is demonstrated here for cancer cell types 
having heightened sensitivity to death receptor activation but also high Bcl2 to Bax ratio (-25). 
Overexpressions for the following signaling molecules are assumed: Bid (2 fold), Bax (2 fold), Bcl2 
(50 fold) and XI AP (5 fold). We performed in silico experiments for the following chemo therapeutic 
targeting of cancer cells: death ligand induction and various degree of Bcl2 and/or XIAP inhibition. 
Increased sensitivity to death receptor activation (upon death ligand induction) was captured by our 
simplified model of DISC clustering (P on = 1, P 0 ff = 10~ 3 , Edd = ~2 K B T). 

DISC formation resulted in significant caspase 8 activation (-10 nM). However, XIAP over-expression 
led to caspase 3 inhibition and strong suppression of apoptotic activation. We observed some capsase 3 

Q 

activation at long times (5><10 MC steps) that can be characterized as mixed type 1 and type 2 
activation. Inhibiting Bcl2 (50-fold over-expression reduced to 15 -fold or 5 -fold) did not alter the 
strong suppression in apoptotic activation (Figure 16). Combined Bcl2 (50-fold to 15-fold) and XIAP 
inhibition (5 -fold to 3 -fold) resulted in significantly increased apoptotic activation involving both type 
1 and type 2 pathways (Figures lie and 16). However, when XIAP was strongly inhibited, type 1 
activation could proceed in a fast deterministic manner (Figures 12d and 16). Such a strategy might also 
allow to selectively target cancer stem cells having heightened sensitivity to death receptor activation. 

Heightened sensitivity to death receptor activation in cancer cells can be accompanied by 
simultaneous increase in anti-apoptotic factors resulting in a dichotomous membrane proximal 
signaling module. Therefore, in spite of the heightened sensitivity to death receptor activation in cancer 
cells (and possibly also in cancer stem-like cells), inhibitory effects of various anti-apoptotic factors 
may only allow low capsase 8 activation leading to type 2 activation (Figures 11a and 12a). The 
inhibitory effects can be potentiated by membrane -proximal factors, such as, increased expression of 
cFLIP [63], increased amount of decoy receptors [64] or altered partitioning of death receptors in 
membrane domains. However, death ligand induced weak apoptotic activation through the type 2 
pathway can be converted to rapid deterministic activation through the type 1 pathway by eliminating 
the increase in anti-apoptotic factors (in the membrane proximal signaling module) and/or enhancing 



Cells 2013, 2 



383 



death receptor density [5]. It has been shown that cFLIP inhibition can render cancer cells susceptible 
to death receptor mediated apoptotic death [63]. Expression of death receptors can be increased by 
inducing transcription [52], increasing transport to membranes, or even by freeing them from other 
binding partners [65,66]. If XIAP is also over-expressed, then simultaneous activation of death 
receptors and inhibition of XIAP will lead to rapid activation of the type 1 pathway (Figure 16). In normal 
cells lack of sensitivity to death receptor activation will lead to protection. Slight over-expression of 
death receptors in some normal cells might initiate type 2 activation, but slow activation with large 
inherent cell-to-cell variability will eventually protect those cells. Recent study in cholangiocarcinoma 
cells demonstrated that a type 2 to type 1 transition in apoptotic activation could be induced in TRAIL 
resistant cancer cells by inhibiting the hedgehog pathway (hedgehog inhibition simultaneously increase 
the death receptor DR4 and decrease XIAP) [18]. We performed further simulations to verify that a 
combined death receptor activation and XIAP inhibition based strategy works better for cancer cells 
having heightened sensitivity to death receptor activation but also high Bcl-2 to Bax ratio. In such a 
scenario, even when Bcl2 (50-fold over-expression reduced to 15-fold) and XIAP (5 fold over-expression 
to normal) were significantly inhibited, 5-fold over-expression of Bid and Bax was not sufficient to 
activate the mitochondrial pathway (through a Bid-Bax type reaction) [12]. However, depending on the 
cancer cell type as well as pro- and anti-apoptotic protein levels, it needs to be determined whether to 
induce a type 1 or a type 2 [12] or a mixed type 1-type 2 activation. Monte Carlo simulations can be 
carried out to simulate various possible options for inducing apoptotic death in cancer cells and 
determine the optimal strategy. 

Figure 16. Activated caspase 3 (fraction of half-maximal) for four different chemotherapeutic 
strategies involving DISC formation and one of the following: (i) Bcl2 inhibition (50-fold 
to 15-fold), (ii) stronger Bcl2 inhibition (50-fold to 5-fold), (iii) Bcl2 (50-fold to 15-fold) 
and XIAP inhibition (5-fold to 3-fold) and (iv) strong XIAP inhibition (5-fold to normal) 
Cancer cells are assumed to have high over-expression levels for both Bcl2 (~50-fold) and 
XIAP (--5-fold). The following parameters are utilized for the DISC clustering: Pon = 1, 
Poff = 10~ 3 , Edd = ~2 KbT. Data analyzed from 64 single cells (MC runs) are shown for 
three different time instants: (a) 2.5 x 10 7 , (b) 5 x 10 7 and (c) 1 x 10 8 MC steps. 
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4. Conclusions 

In this work, we obtain quantitative information regarding the systems level regulation of the 
type 1/type 2 choice in apoptotic death. Activation level of caspase 8 (initiator caspase) emerges as a 
key regulator of the type 1/type 2 choice and stochastic to deterministic transition in apoptotic 
activation. Our results indicate that increased amount of DISC clustering (presumably by higher 
expression of death receptors and/or lower expression of decoy receptors) results in increased caspase 8 
activation (in a rapid manner). Thus cell types with high death receptor expression are expected to 
activate the type 1 pathway more frequently making them type 1 prone (for that specific death 
receptor) [5]. Conversely, low death receptor expression would make cells more prone to type 2 
activation. However, over-expression of certain molecules and/or inhibition of some others may allow 
activation of the type 1 pathway in type 2 prone cells (or the type 2 pathway in type 1 prone cells). 
We also elucidate various mechanisms for mixed type 1 and type 2 activation. Expression levels of 
various downstream (of caspase 8) signaling molecules and kinetic reaction rates were shown to 
determine the time needed for completion of caspase 3 activation (through the two pathways) and have 
impact on the type 1/type 2 choice. Hence, the type 1/type 2 choice in apoptotic activation is regulated 
at a systems level. Interestingly, our results indicate that the type 1/type 2 choice is linked to cell-to-cell 
stochastic variability in apoptosis activation. We elucidate the following contributions to cell-to-cell 
stochastic variability: (i) time-to-caspase 8 generation (membrane proximal signaling module), 
(ii) time-to-cytochrome c release (pre -mitochondrial signaling module), and (iii) time-to-apoptosome 
formation (post-mitochondrial signaling module). The magnitude of cell-to-cell variability depends on 
the cellular context as well as the mechanism of activation. 

In cancer therapy, it is important to perturb the apoptotic signaling network in such a manner that 
cancer cells are selectively targeted and cell-to-cell stochastic variability in apoptotic activation is 
minimized. Interestingly, these two aspects are intricately linked in apoptosis activation of cancer cells. 
In this study, we show that death receptor activation (such as by death ligands), along with XIAP 
inhibition, can lead to rapid activation of the type 1 pathway selectively in certain cancer cells. Thus, it 
provides a mechanism for type 2/mixed type 1-type 2 (slow and often stochastic) to type 1 (rapid and 
deterministic) transition in cancer cells. Such a chemotherapeutic strategy might work for resistant 
cancer cells equipped with high Bcl2 to Bax ratio but also having increased sensitivity to death receptor 
activation. Normal cells having lower death receptor expression, in contrast, would remain protected 
due to slow activation and inherent variability of the type 2 pathway (under the application of death 
ligand and XIAP inhibitor). Our studies indicate how in silico experiments can be performed to 
determine the optimal strategy for targeting the apoptotic pathway of cancer cells in a cell-type 
specific manner. 

Appendix I 



Table Al. Reaction moves considered in our Monte Carlo simulation. 



Apoptosis signaling reactions 


Kinetic rate constants 




DISC + ProCaspase8 ;± Disc-ProCaspase8 


kon: 


3.5 x lO^nlvr 1 s _1 ;k off : li 


S x 10 2 s 1 [6] 


Disc-ProCaspase8 + ProCaspase8 ^ Disc-ProCaspase8_2 


kon- 


3.5 x KT'nlvr 1 s _1 ;k off : U 


i x 10 2 s 1 [6] 


Disc-ProCaspase8 2 DISC + 2 P 4 3/ 41 


k C at- 


3.0 x 10~' s" 1 [6] 
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Table Al. Cont. 



Apoptosis signaling reactions 

P43/41 -» Caspase 8 

Caspase 8 + Bid ^ Caspase8-Bid 

Caspase8-Bid -> Caspase 8 + tBid 

Caspase 8 + ProCaspase 3 ^ Caspase8-ProCaspase 3 

Caspase8-ProCaspase 3 -> Caspase 8 + Caspase 3 

Bid + Bax £ Bid-Bax 

Bid-Bax + Bax «» Bid-Bax 2 

Bid-Bax 2 -> Bid + Bax 2 

Bid + Bcl 2 2 Bid-Bcl 2 

tBid + Bax ♦» tBid-Bax 

tBid-Bax + Bax £ tBid-Bax 2 

tBid-Bax 2 -> tBid + Bax 2 

Bcl 2 + tBid *» Bcl 2 -tBid 

Bcl 2 + Bax ♦» Bcl 2 -Bax 

Cytochrome C + Apaf *» Cytochrome C-Apaf 

Apaf + Apaf 5± Apaf- Apaf 

XIAP + Smac <± XIAP-Smac 

Apoptosome + ProCaspase9 *i Apoptosome-ProCaspase 9 
Apoptosome-ProCaspase9 + ProCaspase9 ^ 
Apoptosome-ProCaspase 9 2 

Apoptosome-ProCaspase 9 2 -> Apoptosome + 2 Caspase 9 

XIAP + ProCaspase 9 ♦» XIAP-ProCaspase 9 

XIAP + Caspase 9 ♦» XIAP-Caspase 9 

Caspase 9 + ProCaspase 3 5^ Caspase 9-ProCaspase 3 

Caspase 9-ProCaspase 3 -> Caspase 9 + Caspase 3 

XIAP + Caspase 3 g XIAP-Caspase 3 



Kinetic rate constants 

k cat : 1.0 x 10 1 s 1 [6] 

k on : 5.0 x 10~ 3 nlVT 1 s" 1 ; k off : 5.0 x 10~ 3 s" 1 [6] 
k cat : l.OxlO'V 1 [6] 

k on : 1.0 x 10~ 5 nlVT 1 s" 1 ; k off : 6.0 x 10~ 2 s~ ! [6] 
k cat : 1.0 x 10 1 s 1 [6] 

k on : 1.0 x lO^nftT 1 s" 1 ; k off : 1.0 x 10~ 2 s~ ! [12] 
k on : 1.0 x lO^nlVT 1 s" 1 ; k off : 1.0 x 10~ 2 s" 1 [12] 
k cat : 1.0xl0" 2 s _1 [6] 

k on : 1.0 x lO^nM -1 s" 1 ; k off : 1.0 x 10~ 2 s" 1 [12] 

k on : 2.0 x lO^nftT 1 s" 1 ; k off : 2.0 x 10~ 2 s" 1 [6] 

k on : 2.0 x lO^nNT 1 s" 1 ; k off : 2.0 x 10~ 2 s~ ! [6] 

k cat : 2.0 x 10 2 s 1 (taken as ~k of f) 

k on : 2.0 x 10~ 3 nlVT 1 s" 1 ; k off : 2.0 x 10~ 2 s" 1 [6,67] 

k on : 2.0 x 10~ 3 nlVT 1 s" 1 ; k off : 2.0 x 10~ 2 s" 1 [6,67] 

k on : 2.8 x 10~ 7 nlVT 1 s" 1 ; k off : 5.7 x 10~ 3 s _1 [6] 

k on : 1.0 x 10~ 2 nlVT 1 s" 1 ; k off : 10 s" 1 

(|xM affinity based on CARD-CARD interaction [7]) 

k on : 7.0 x 10~ 3 nlVT 1 s" 1 ; k off : 2.2 x 10~ 3 s" 1 [6,68] 

k on : 2.8 x lO^nlVT 1 s" 1 ; k off : 7.5 x 10~ 2 s~ ! [6] 

k on : 2.8 x lO^nftT 1 s" 1 ; k off : 7.5 x 10~ 2 s _1 [6] 

k cat : 7.0 x I0 l s" 1 [6] 

k on : 5.0 x 10~ 3 nlVT 1 s" 1 ; k off : 3.5 x 10~ 3 s" 1 

(taken to be same as that for Caspase 9) 

k on : 5.0 x 10~ 3 nlVT 1 s" 1 ; k off : 3.5 x 10~ 3 s _1 [25,69] 

k on : 2.0 x 10~ 5 nNT 1 s" 1 [6,26]; k off : 5.7 x 10 _1 s" 1 [6,25] 

k cat : 4.8 s" 1 [6,70] 

k on : 2.5 x 1(T 3 nNT 1 s" 1 ; k off : 2.4 x 1Q' 3 s ~' [6,71] 



Table A2. Initial number of molecules. 



Molecule type Concentration (nM) 

DISC (adaptor such as FADD/TRADD) 1 00 [32] 
ProCaspase 8 30 [3,6,72] 

Bid 33 [6,25] 

Bax 83 [6] 

Bcl-2 75 [6] 

Cytochrome C 100 [6] 

Smac 50 

(-100 nM in [6]) 

Apaf 100 [6]; 20 [26] 
XIAP 30 [3,6] 
ProCaspase 9 20 [6,70,73] 
ProCaspase 3 100 [72] 
(also see [3,6,23]) 
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Table A3. Probabilistic simulation parameters. 



Apoptosis signaling reactions 



Probabilistic parameters used in 
the Monte Carlo simulation of 
apoptosis signaling (P on /Pofl/Pcat) 



DISC + ProCaspase8 □ Disc-ProCaspase8 

Disc-ProCaspase8 + ProCaspase8 □ Disc-ProCaspase8 2 

Disc-ProCaspase8_2 -> DISC + 2 P43/41 

P43/41 -> Caspase 8 

Caspase 8 + Bid □ Caspase8-Bid 

Caspase8-Bid — > Caspase 8 + tBid 

Caspase 8 + ProCaspase 3 □ Caspase8-ProCaspase 3 

Caspase8-ProCaspase 3 — > Caspase 8 + Caspase 3 

Bid + Bax □ Bid-Bax 

Bid-Bax + Bax □ Bid-Bax 2 

Bid-Bax 2 -> Bid + Bax 2 

Bid + Bcl 2 □ Bid-Bcl 2 

tBid + Bax □ tBid-Bax 

tBid-Bax + Bax □ tBid-Bax 2 

tBid-Bax 2 -» tBid + Bax 2 

Bcl 2 + tBid □ Bcl 2 -tBid 

Bcl 2 + Bax □ Bcl 2 -Bax 

Cytochrome C + Apaf □ Cytochrome C-Apaf 
Apaf + Apaf □ Apaf-Apaf 
XIAP + Smac □ XIAP-Smac 

Apoptosome + ProCaspase9 □ Apoptosome-ProCaspase 9 

Apoptosome-ProCaspase9 + ProCaspase9 □ Apoptosome-ProCaspase 9_2 

Apoptosome-ProCaspase 9 2 -> Apoptosome + 2 Caspase 9 

XIAP + ProCaspase 9 □ XIAP-ProCaspase 9 

XIAP + Caspase 9 □ XIAP-Caspase 9 

Caspase 9 + ProCaspase 3 □ Caspase 9-ProCaspase 3 

Caspase 9-ProCaspase 3 -> Caspase 9 + Caspase 3 

XIAP + Caspase 3 □ XIAP-Caspase 3 



P on : 1.0; P off : 1.8 x 10~ 6 

(P on taken as 1.0 instead of 0.35 as 

an approximation) 





1.0; P 


off : 1.8 x 10~ 6 






Peat 


3.0 


X 


10" 5 








Peat 


1.0 


X 


10- 5 








P * 
A on- 


5.0 


X 


lO^Poff 


5.0 


x 10 


-7 


Peat 


1.0 


X 


10- 5 








Pon- 


1.0 


X 


10~ 3 ; Poff 


6.0 


x 10 


-6 


Peat 


1.0 


X 


10" 5 








P • 

A on* 


1.0 


X 


10~ 4 ; P off : 1.0 


x 10 


-6 


p • 

A on* 


1.0 


X 


10~ 4 ; P off : 1.0 


x 10 


-6 


Peat 


1.0 


X 


10" 6 








Pon* 


1.0 


X 


10~ 3 ; Poff 


1.0 


x 10 


-6 


Pon* 


2.0 


X 


10~ 2 ; Poff 


2.0 


x 10 


-6 


Pon- 


2.0 


X 


10~ 4 ; Poff 


2.0 


x 10 


-6 


Peat 


2.0 


X 


10~ 6 








A on- 


2.0 


X 


10 _1 ;P off 


2.0 


x 10 


-6 


Pon- 


2.0 


X 


10 _1 ;P off 


2.0 


x 10 


-6 


Pon* 


2.8 


X 


10~ 5 ; Poff 


5.7 


x 10 


-7 


Pon* 


1.0 


X 


10°; P off : 1.0 x 10~ 3 


Pon* 


7.0 


X 


lO'^Poff 


2.2 x 10 


-7 


Pon- 


2.8 


X 


10" 2 ; Poff 


7.5 


x 10 


-6 


P * 
A on- 


2.8 


X 


10~ 2 ; Poff 


7.5 


x 10 


-6 


Peat 


7.0 


X 


10 5 








Pon- 


5.0 


X 


10 _1 ;P off 


3.5 


x 10 


-7 


Pon* 


5.0 


X 


10"'; Poff 


3.5 


x 10 


-7 


Pon* 


2.0 


X 


10~ 3 ; Poff 


5.7 x 10 


-5 


Peat 


4.8 


X 


10- 4 








P * 

A on- 


2.5 


X 


10 _1 ;Poff 


2.4 x 10 


-7 



Appendix II 

The kinetic rate equations for active caspase 8 reacting with procapsase 3 (type 1 pathway) and Bid 
(type 2 pathway) are provided below [6,23,25]. [] denotes concentration of a given signaling molecule. 
We use the following abbreviations: C8 (active caspase 8), C3 (pro caspase 3), C3 * (active caspase 3). 
d[C8\ 

= -fc^ n [C8][C3] + kl ff [C8_C3] + k^ t [C8_C3] - [C8][Bid] + k\ n [C8_Bid] + k^ t [C8_Bid] (1) 
dlC3J 

(2) 
(3) 
(4) 



d[C8-C3l 
dt 



= -4, [C8][C3] + kl f f [C8_C3] 

= k^ [C8][C3] - k$f f [C8_C3] - kt at [C8_C3] 
d[C3*\ 



dt 



= kLr[C8_C3] 



Cells 2013, 2 



387 



dt 

d[C8-Bid\ 



= -fc*. [C8] [Bid] + kl Jf [C8_Bid] (5) 



= *=, [C8] [Bid] - kl ff [C8_Bid] - kl ct [C8_Bid] (6) 
= khtlCQ.Bid] (7) 



An approximate estimation of caspase 8-procaspase 3 and caspase 8-Bid complexes is carried out 
by (i) assuming that the system has reached a steady state and (ii) ignoring the catalytic activation 
terms (k cat y. 

[C8.C3] =k^[C8][C3] (8) 
[C8_Bid] = ki[C8][Bid] W 

where Ka = k on lk 0 ff 

[C8_Bzd] _ k^[C8][Bid] _ k^iBid 0 - [C8_Bid]) 

~CS_C3] ~ k-:C8Y£3} ~ - ;C8_C3] ( 10 ) 

C3° and Bid 0 denote the initial concentrations of procaspase 3 and Bid, respectively. For low caspase 
8 concentrations (~1 nM), [C8 C3] « C3° (100 nM) and [C8_Bid] « Bid 0 (33 nM). In such a limit: 

[C8 sii ] kjj[C8][Bid] kjBid* 

[C8„] " fc^[C8][C3] " k* C3 S " 10 ' ( U ) 
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